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Transport mechanism in the normal state of high-T c superconductors is described using the well 
known Fermi-Dirac statistics in which an additional restrictive constraint is introduced so as to 
capture the variation of resistivity with temperature and doping. The additional restrictive condi- 
tion is the ionization energy that will eventually determine the properties of charge carriers' in the 
normal state of high-T c superconductors. The magnitude and the variation of charge carriers con- 
centration and resistivities (polycrystalline, c-axis and afe-planes) with temperature and doping are 
very well described by the ionization energy based Fermi-Dirac statistics. However, these transport 
models are not appropriate for cuprates below the characteristics (T*) and critical temperatures 
(T c ), metals with free electrons and strong electron- phonon scattering. Ionization energy is found 
to be an essential parameter to accurately predict variations of resistivity's magnitude with doping, 
charge carriers' concentrations, scattering rate constants as well as the effective mass. Apart from 
that, iFDS based resistivity models provide the comprehensive information on crossover-temperature 
(metallic — > insulating transition temperature) in the normal state of high-T c superconductors. 

PACS numbers: 71.10.Ay, 74.72.Bk, 72.60.+g, 74.72.-h 



I. INTRODUCTION 



Superconductivity has been around for almost a century since its discovery in 1911 by Kamerlingh-Onnes [1]. 
Subsequent discovery on Copper-oxide (CU-O2) materials by Bednorz and Muller [2] in 1986 literally questions the 
applicability of the Nobel prize winning Bardeen-Cooper-Schrieffer (BCS) theory [3] in cuprates. This terrible yet 
exciting scenario have had led to numerous theoretical proposals and experimental endeavours due to the fact that the 
transport properties of CU-O2 based ceramics are overwhelmingly mystifying. Apart from cuprates, superconductivity 
seems to occur in almost all type of materials namely, organic Bechgaard salts [4, 5], ferromagnetic [6] ZrZn2, metallic 
Magnesium diboride [7], Buckminstcrfullcrcncs [8], Carbon- nanotubes [9] and unbelievably even in DNA [10]. It seems 
that superconductivity could be one of the fundamental properties of nature that has been elusive prior to Bednorz's 
discovery. Cuprates and its potential applications were also extensively publicized to the extent where it initiated the 
superconductivity-renaissance, which was between late 80s and early 90s due to exceptional increase of T c from about 
20 K to 90 K within months. In those days before 1986, reports on superconductivity were revolved around simple 
metals, metallic alloys and some organic salts. Perhaps due to bad omen, say brought about by the black schrodinger 
cat, scientists have had a hard time to find the correct mechanism for conventional superconductors, which nearly took 
them 50 years to arrive at BCS theory in 1957. Hopefully history will not repeat itself for the superconducting cuprates 
that further perplex scientists with its puzzling phenomena. Cuprates, the unique oxide has many faces depending 
upon temperature and doping unlike other known superconductors and electronic materials reported thus far. One 
of the face is of course, superconducting region but there are more than meets the eye such as antiferromagnetic 
fluctuations, pseudogap, Fermi liquid and non-Fermi liquid as one dopes the cuprates or varies the temperature. As 
a matter of fact, there have been sufficient amount of experimental results on cuprates that could somehow pin-point 
the proper ingredients for at least to understand the normal state transport properties of high-T c superconductors 
to a certain reasonable extent. Anisotropic phenomena on the normal state resistivity (p(T)) curves such as the 
existence of metallic-like conduction, semiconducting property below a certain crossover-temperature {T crossover ) and 
incompatibility of T crossover in c-axis with characteristic temperature (T*) in afc-plane [11, 12] point towards that 
Tcrossover in Pc(T) could be a different entity from T* in p a b{T). Therefore, it is possible to work-out the details of 
transport properties in c-axis in accordance with ionization energy based Fermi-Dirac statistics (iFDS) in which the 
charge carriers is further confined with an additional parameter known as the ionization energy (Ej). In the procedure 
of deriving iFDS, a somewhat different distribution functions for holes and electrons are also obtained that further 
justifies the inappropriateness of iFDS for systems consist of free electrons and also for systems which are influenced 
by strong electron-phonon scattering or both [13, 14]. It is a well established fact that the transport dimensionality in 
the normal state of high-T c superconductors is 2-dimensional (2D). As such, comprehensive derivation of iFDS will be 
carried-out to model the 2D equations of charge carriers' concentrations {y/rvp(T,Ej) and n(T,Ej)) and resistivities 
{PpoiyiT, Ej) and p c (T, Ei)) as a function of temperature (T) in the vicinity of electric field that will eventually enable 
one to enumerate y/np(T, Ej), n(T, Ej), p po i y {T, Ej) and p c (T, Ej). Gapless nature of charge carriers in afe-planes as 
a consequence of the experimental results by Basov et al. [15] is due to an interesting discrepancy between a6-plane 
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and c-axis in term of supcrfluid density (p s ) and spectral weights in the normal (N n ) and superconducting states 
(N 3 ). In c-axis, p s > N n — N s while p s = N n — N s in afr-plane, which signifies a gaplcss behavior of charge carriers 
in afe-planes [16]. This result will also be incorporated into p po i y (T,Ej) to obtain p a b{T). Previously, Anderson and 
Zou [17] proposed an empirical resistivity equation given by 



P = ^ + BT, (1) 

for both a6-plane and c-axis. Equation (1) simply suggests that there are two types of competing T-dcpcndcnces in 
the normal state of high-T c superconductors. It is quite obvious that the coexistence of metallic and semiconducting 
resistivities in both p ao (T) and p c (T) are intriguing. Basically, p c (T,Ei) transition occurring from metallic (normal 
state) superconductor (MS) to metallic-semiconducting crossover (normal state; without T*) superconductor (MSS) or 
vice versa with doping is shown to be predictable with models derived from iFDS. It is to be noted that the significance 
of T* on p a b(T) is ignored in this work since it is still controversial both in terms of experimental and theoretical 
approaches [11, 12]. Parallel to this, p a b(T) and R^ b \T) involving T* are yet another critical sub-aspects of the 
above-mentioned conduction peculiarities that need further studies. A theoretical proposal on these phenomena that 
consist of electrons, spinons and holons was given in Rcfs. [11, 12] and references therein. Finally, superconducting 
samples of both crystalline and polycrystalline reported by various research groups will be employed extensively to 
accentuate the importance of iFDS in the normal state of 2D high-T c superconductors. 



II. IONIZATION ENERGY BASED FERMI-DIRAC STATISTICS 



The conduction electron's distribution can be derived using iFDS with ionization energy as an anomalous constraint. 
This derivation involves two restrictive conditions: (i) the total number of electrons in a given system is constant and 
(ii) the total energy of n electrons in that system is also constant. Both conditions are as given below 



^2 n t = n, (2) 



^TEim = E. (3) 

i 

However, the condition as given in Eq. (3) must be rewritten as given in Eq. (4) by inserting conditions, E e i ectron 
- Einitiaistate + Ei and Ehoie = E init i a i state — Ei appropriately. 



y^^(Ei ni tialstate ± E^jTlj — E. (4) 

i 

This is to justify that an electron to occupy a higher state N from initial state M is more probable than from 
initial state L if condition Ej(M) < Ej(L) at certain T is satisfied. As for a hole to occupy a lower state M from 
initial state N is more probable than to occupy state L if the same condition above is satisfied. E ini tiai state {Ei. s ) is 
the energy of a particle in a given system at a certain initial state and ranges from +00 to for electrons and to 
—00 for holes. The importance of this inclusion is that it can be interpreted as a gap that will be described later and 
also, particularly the Ej can be used to estimate the resistivity transition upon substitution of different valence state 
ions. Utilizing Eq. (4), one can write the mean number of particles in a certain quantum state s (n s ) as 

E n, exp[-A[ni(£7 ia ± £1)1 + n 2 (E is ± Ej) 2 + ... + n s (E is ± Ei) s + ...]] 
E exp[-A[m (E is ±E I ) 1 +n 2 {E is ±E I ) 2 + ... + n s (E is ± Ej) s + ...]} y > 

Since state s was chosen arbitrarily to calculate the mean number of particles in it thus Eq. (5) has to be algebraically 
rearranged as given below to ease the following derivation 
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^n s exp[-\n s (E is ±E I ) s ] £ exp[-A[m(-E is ± + n 2 (E is ± £/) 2 + ...] 

n 3 ni,n 2 ,... 

J2exp[-Xn s {E is ±E I ) s ] £ exp[-A[n 1 (^ s ±S / ) 1 +n 2 (S JS ±S / ) 2 + ...]] 



n s ni,n2.... 

Note that the sum, X^n 2 omits the chosen state, s. Subsequently, the partition function that represents 
Sni ri2 can ^ e wr itten as 



Z a (n)= J! exp[-A[ni(^ s ± J B 7 ) 1 + n 2 (E is ±E I ) 2 + ...}]. (7) 



Fermi-Dirac statistics, including iFDS also need to satisfy Pauli's exclusion principle apart from the restrictive 
conditions given in Eqs. (2) and (4) in which, Pauli's principle requires that one needs to sum over n s = and 1. This 
implies that each state can only accomodate one particle. In doing so, one arrives at Eq. (8) from Eqs. (6) and (7). 

_ + cxp[-\n s (E is ±E I ) s ]Z s (n-l) 



Z s {n) + exp[-Xn s (E is ± E I ) s ]Z s (n - 1) 

(8) 



^Lexp[An s (£; is ±£ 7 ) s ] + r 
Further simplification of Eq. (8) can be performed using the relation (assuming n»l) 

\nZ s {n-l) = \aZ s {n)- ° Zf 



Z s (n — 1) = Z s (n) exp 



dn 
d In Z s 



<9n 

Taking din Z s /dn = fi s =(iasa consequence of large n and substituting it into Eq. (8) will lead one to write 



(9) 



1 

ns " e X p\ M + X(E ia ±E I ) s ] + l' (10) 

Equation (10) is the ionization energy based Fermi-Dirac statistics and the procedure for deriving it, is identical to 
standard Fermi-Dirac statistics that is given in Ref. [18]. 

III. RESISTIVITY MODELS 

The probability functions for electrons and holes are derived from Eq. (10) after further considering exp[^j,+X(E±Ei)] 
^> 1 which are given in Eqs. (11) and (12) respectively [16] 

f e {E,E I )=exp[-n-\(E + E I j\, (11) 

fh(E, Ej) = exp [/i + \(E — Ej)j . (12) 

The parameters (x and A are the normalization constants that will be derived shortly from Eqs. (2) and (4) respec- 
tively In the standard FDS, the respective Eqs. (11) and (12) are simply given by 



f e {E) = exp[-fx-\E], 



(13) 
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f h (E)=e X p[n + XE]. (14) 
The 2D density of states' (DOS) derivative can be shown to be 

^ = ^. (15) 

Subsequently, the restrictive condition as given in Eq. (2) can be rewritten for both electrons and holes respectively 
by employing the above-stated 2D dN/dk as 

„= -e-^ykexp^A—jdk, (16) 
o 



,-£-»./ k «p(Ai£)* (17) 

— oo 

Note here that is substituted with h 2 k 2 /2m. m is the mass of the charge carriers while n and p are the respective 
concentrations of electrons and holes. L? denotes area in k-space and h — h/2n, h is the Planck constant. The 
respective solutions of Eqs. (16) and (17) are given below 

- ^ < M » 

The above-stated solutions were obtained via integral 

oo 

/ x 2 " +1 cxp(-ax 2 )dx= -^L. 
o 

Equations (18) and (19) respectively imply that [14] 

He(iFDS)= He + \E I: (20) 

l i h {iFDS)=n h -XE I . (21) 
On the other hand, using Eq. (4), one can obtain 

2fi2 °? / , fi 2, 2 - 



£ = e -**-^ / k 3 cxp rfk 



m f L 
2tt V Aft" 



o 

2 

e^- AE '. (22) 



Equation (22) after appropriate substitution with Eq. (18) is compared with the energy of a 2D ideal gas given 
below 
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E = nk B T. 



(23) 



Quantitative comparison will eventually leads to \{fds = ^fds = 1/ksT where k B is the Boltzmann constant. 
Actually, equation (23) can be derived trivially with the procedure given in Refs. [18, 19]. Identical A between FDS and 
iFDS has been somewhat anticipated and is not unusual cither. Nevertheless, Eqs. (20) and (21) will have enormous 
consequences on the properties of charge carriers as will be pointed out in the following paragraphs. 

The distribution function for electrons and holes can be written explicitly by first denoting fj, = —Ep (Fermi level), 
A = l/k B T and substituting these into Eqs. (11) and (12) will lead one to write 



f e (E,Ei) = exp 



Ep — Ej — E 



(24) 



fh(E,E T ) = exp 



E — Ej — Ep 



knT 



(25) 



It is worth noting that, — Ej in Eq. (25) for holes is a natural consequence of Dirac's theory of antiparticle inter- 
pretations [20] . Physically, Ej is due to electron-ion attraction or Coulombic in nature that has to be contrasted with 
band gap, E g which arises from energy band splitting. General equations to compute charge carriers' concentrations 
are as listed below 



oc 

n = J f e (E,E I )N e (E)dE. 
o 



(26) 



u 

V = J fh(E, E I )N h (E)dE. 



(27) 



N e (E) and Nh(E) are the 2D DOS for electrons and holes respectively. It is now convenient to obtain the geometric- 
mean concentrations of electrons and holes in the normal state of 2D polycrystallinc high-T c superconductors as [16] 



np(T,Ei) = -^(m* e m* h ) 1/2 exp 
irh 

Here, n w p is assumed to avoid Ep and the 2D DOS is given by 



-Ej 

k B T 



(28) 



N e , h (2D) 



"e,h 



irh 



2 ■ 



(29) 



Note here that iV ej /j(2D) is independent of E and m* h is the effective mass of either electrons or holes. In order to 
calculate p c (T 7 Ej), n needs to be determined first that is given by (from Eq. (26)) 



n(T, Ej) 



m*kBT 
5 — exp 

TTh 2 



Ep — Ei 



k B T 



(30) 



The 2D resistivity models for polycrystals and in c-axis can be derived from an elementary resistivity equation, 



(31) 



Substituting, 1/r = AT 2 as a consequence of electron-electron interactions and Eq. (28) for n into Eq. (31) gives [16] 
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Ppoly 

Similarly, p c (T 7 Ej) can be shown to be [13] 



e 2 k B 



Texp 



111 

k B T 



(32) 



(33) 

Note that A 2 is a 2D T-independent scattering rate constant in the absence of magnetic field and e is the charge of 
an electron. Furthermore, one can recall the characteristics of Eq. (1) which has been captured by Eqs. (32) and (33) 
where both equations also contain two competing T-dcpendcnces though in a different proportionality. One can also 
surmise that for 2D systems, k B T crossover = E I for p po i y (T) and k B T crossover = Ei - E F for p c {T). Now, applying 
the result of gapless nature of charge carriers in a6-planes into Eq. (32) leaves one with [16] 



Pc(T,Ej) 



A 2 irh z , 
e 2 k B 



■Texp 



Ej — Ef 
k B T 



Pab(T) 



e 2 k B 



T. 



(34) 



The nature of gapless phenomena in a6-planes could be due to change in scattering rate that originates from a 
different Fermi phase space argument in accordance with nested Fermi liquid theory [21]. Note that Eq. (34) is 
nothing but a special case that satisfy, Ej <C T c inequality It is stressed that the influence of free electrons and 
T-dependence of electron-phonon scattering were entirely ignored in the p po i y (T, Ej) and p c (T, E[) models as given in 
Eqs. (32) and (33) respectively. Hence, those models are obviously not suitable for conventional metals. Interestingly, 
the electrons of a 2D semiconductor with E g can be further gapped with Ei that will lead one to derive v / np(T, Ej,E g ) 
and p(T, Ei, Eg) respectively as 



n P(T, Ei,E g ) 



k B T 
nh 2 ' 



* *\!/ 2 

m e m h ) exp 



-Ei 



2 E 9 



k B T 



(35) 



p(T,Ei,E g ) 



A 2 nfr 
e 2 k B 



Texp 



Ei 



2 E 9 



k B T 



(36) 



Consequently, electrons in the conduction band (E > E g ) will still be influenced by Ei even though Ei < E g due 
to Ei's Coulombic nature between charge carriers and ions. 



IV. DISCUSSION 



The effect of Nd 3+ (Tj = 1234 kJmor 1 ) substitution into Sr 2 + (Ei = 807 kJmor 1 ) in superconducting 
TlSr2_ 2: Nd a ;CaCu207 compound [22] was found to increase the p(T) in accordance with Ei. This justifies the 
need for Ei based analysis on doping as pointed out by iFDS. Applications of iFDS in superconductors are ex- 
plicitly given in Rcfs. [11, 12, 13, 16]. Recently, Naqib et al. [23] have investigated the electrical properties of 
Yi_ a; Ca 2 ;Ba2(Cui_j / Znj ) )307_ ( j superconducting compounds by varying x, y and d. The transition of normal state 
p(T) with Ca 2+ (Ei = 867 kJmol -1 ) doped into Y 3+ (Ej = 1253 kJmol -1 ) sites is in excellent agreement with 
Eq. (32) of iFDS. But, Zn 2+ doping is not appropriate to analyze as a function of iFDS only because this substi- 
tution will directly disturb the a6-plane conduction of spinons and holons and also in term of oxygen concentration 
(d), thus the overall conductivity of Yi_ a; Ca a ;Ba2(Cui_ y Znj / )307_ ( i polycrystals will be modified in a not-so-simple 
way [11, 12]. It is easy however, to extract the relation of normal state p(T) between Yo.9Cao.iBa2Cu307_ ( z and 
Yo.8Cao.2Ba2Cu3C , 7_ < i where the normal state p(T) is reduced with Ca 2+ doping for all d (oxygen pressure), since 
Y 3 +(Ei = 1253 kJmol- 1 ) > C& 2 +(Ei = 867 kJmor 1 ). 

In contrast, Sr 2+ (Ei = 807 kJmol -1 ) substitution into Ba 2+ (Ei = 734 kJmol -1 ) sites have decreased the normal 
state p(T) in Hg .85Reo.i5(Bai_ y Sr y )2Ca2Cu308-,5 unexpectedly [24]. iFDS suggests that p(T) should increase with 
Sr doping into Ba sites since Sr 2+ (Ei = 807 kJmol -1 ) > Ba 2+ (Ei = 734 kJmol -1 ). This contradicting scenario can 
be explained since the actual doping concentration determined with EDX showed that the concentrations of other 
elements namely, Re, Ca and Cu also vary with Sr doping into Ba sites. It is quite non-trivial to verify and prove 
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these by calculating the relative Ej for Ca 2+ , Sr 2+ , Ba 2+ , Re and Cu 2+,3+ . For comparison purposes, the valence 
state of Re and Cu are taken to be 3+ and 2+ respectively. Arbitrary values of valence state for Re and Cu will 
not affect this analysis since its valence states are assumed to be constant (in this case) for SrOO, Srl7 and Sr28. It 
is known from Ref. [24] that the concentrations and Ej for SrOO are Re 3+ (0.15; 1510 kJmor 1 ), Ba 2+ (2.10; 734 
kJmor 1 ), Sr 2 + (0.00; 807 kJmor 1 ), Ca 2+ (2.20; 867 kJmor 1 ) and Cu 2 + (3.10; 1352 kJmor 1 ). For Srl7, it is given 
by Re 3+ (0.15; 1510 kJmor 1 ), Ba 2 + (0.84; 734 kJmor 1 ), Sr 2 + (0.17; 807 kJmor 1 ), Ca 2 + (1.97; 867 kJmor 1 ) and 
Cu 2 + (3.12; 1352 kJmol- 1 ). Finally for Sr28, Re 3 + (0.14; 1510 kJmol- 1 ), Ba 2 + (0.74; 734 kJmol- 1 ), Sr 2 + (0.28; 
807 kJmor 1 ), Ca 2+ (1.75; 867 kJmor 1 ) and Cu 2+ (3.02; 1352 kJmor 1 ). Therefore, from this data it is possible to 
calculate the changes of Ej due to the fluctuations of other non-dopant elements' concentrations with Sr doping into 
Ba sites. One can show that the relative Ejs are as given below for SrOO, Srl7 and Sr28 respectively. 



E\ b rUUJ = [0.15(1510)] Re3+ + [2.10(734)] Ba 2+ 
+ [0.00(807)] Sr 2+ + [2.20(867)] Ca 2+ 
+ [3.10(1352)] Cu2+ 
= 7640 kJmor 1 , (37) 



E\ brl7) = [0.15(1510)] Re3+ + [0.84(734)] Ba2 + 
+ [0.17(807)] Sr 2 + + [1.97(867)]c«»+ 
+ [3.12(1352)] Cu2+ 
= 6680 kJmor 1 , (38) 



E<f r28) = [0.14(1510)] Re3+ + [0.74(734)] Ba = + 
+ [0.28(807)] Sr 2 + + [1.75(867)] Ca »+ 
+ [3.02(1352)] Cu2+ 
= 6369 kJmor 1 . (39) 

Hence, the reduction of p(T) with Sr doping is justified from Eqs. (37), (38) and (39), which is due to the concen- 
tration's fluctuation of Ca, Re and Cu apart from Ba and Sr. All values of Ej including in Eqs. (37), (38) and (39) 
were averaged in accordance with 



i=i z 

and should not be taken literally since those Ejs are not absolute values. The absolute values need to be obtained 
from Eq. (41) below [14] 

Ej = JT^— (41) 

If one of the substituting ion is multivalcncc, then the valence state of that particular ion can be approximated 
using Eq. (42) below [16] 

- E En + ^En = -±E H . (42) 

J i=z + l i=l y i=l 

The variations of p(T) with doping in polycrystallinc samples namely, (Pri_ x Gd x )i.85Ceo.i5Cu04 and 
(Pri_a;Y x )i.85Ceo.i5Cu04 from Meen et al. [25], Pr :c Gdi_ a ;Ba2Cu307_5 from Khosroabadi et al. [26] can be explained 
in a straight-forward manner without even using Eq. (42). Meen et al. [25] found that substitutions of Gd 3+ and Y 3+ 
into Pr 3+ give rise to p(T) simply because £/(Pr 3+ = 1211 kJmor 1 ) < £/(Gd 3+ = 1251 kJmor 1 ) and £ 7 (Pr 3+ ) 
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< Ej(Y 3+ = 1253 kJmol -1 ). Similarly, Khosroabadi et al. [26] observed decreasing p(T) for Pr 3+ substitution into 
Gd 3+ due to former inequality. Further indirect observation was given by Isawa et al. [27] for Ce ;c Pr 1 _ a ,La2Cu04 +( 5 
polycrystals. However, the valence state of Ce is unknown in that sample therefore one has to employ Eq. (42). The 
first term, | X)i=f+i m Eq- (42) has i = z + 1, z + 2,..., z + j and j = 1, 2, 3,.... It is solely due to multivalence 
ion. In this case Ce 3+ ' 4+ is substituted into Pr 3+ sites hence from Eq. (42), the first term is due to Ce 4+ ion's 
contribution or caused by reaction of the form Ce 3+ — electron — > Ce 4+ (3547 kJmol -1 ), hence j is equals to 1 here 
and S represents the additional contribution from Ce 4+ . The second (i — 1, 2, 3, z) and last (i = 1, 2, 3, q) 
terms respectively are due to reaction of the form Ce — 3 (electrons) — ► Ce 3+ and Pr — 3 (electrons) — ► Pr 3+ . Recall 
that q = z = 3+ and i = 1, 2, 3,... represent the first, second, third, ... ionization energies while j = 1, 2, 3, ... 
represent the fourth, fifth, sixth, ... ionization energies. Therefore, z + S gives the minimum valence number for Ce 
which can be calculated from Eq. (42). Apparently, Ce 3+ [Ej = 1178 kJmol -1 ) substitution into Pr 3+ (Ej = 1211 
kJmol -1 ) will lead to a lower p{T) as indicated in Ref. [27]. The opposite scenario, increment of p(T) with Ce 3+<5 is 
possible if and only if the valence state of Ce is > 3.010. This value (3.010) is actually obtained from Eq. (42). 

Interestingly, Shi et al. [28] have reported substitution effects of Rui-sSb^S^Gdi^Ceo.eCu^Oio-a polycrystals on 
p(T). Since £/(Sb 3+ = 1623 kJmor 1 ) < £/(Ru 3 + = 1692 kJmor 1 ) then it is quite obvious to predict p{T) will be 
reduced in magnitude with Sb 3+ content. However, a contradicting scenario have been observed where p(T) increases 
with Sb. As such, Eq. (42) is again useful to estimate the minimum valence state of Sb 3+<5 in order to justify p(T) 
increment with Sb. The minimum valence state of Sb has to be > 3+ (z + 5) so as to satisfy p(T) increment. 
z + 5 is calculated to be 3.017 or Sb 3 017+ . Another point worth extracting is that additional annealing (4 days) 
for RuSr 2 Gdi.4Ceo.6Cu20io-<5 (sample B) does not change Ej significantly whereas the A 2 parameter is reduced 
dramatically compared to as-prepared sample (A). The assumption that Ej changes insignificantly with annealing 
is valid since the shape of p(T) curves for samples A and B are identical [28]. This phenomenon clearly indicate 
that long hours of proper heat-treatment improves sample B's quality in terms of grain boundaries, defects [28] 
and subsequently reduces the magnitude of scattering rates (A 2 ). Figure 1 indicates calculated plots of p(T,Ej) for 
Hgo.7Cr .3Sr 2 Cu04 +( 5 from Ref. [29] and Gdo.95Pr .o5Ba2Cu307_5 from Ref. [26] while the inset is for Sr .gLao.iCu02 
from Ref. [30] and TmBa2Cu3 6.99 from Ref. [31]. Note that the symbols do not represent the experimental data 
points but simply to enhance contrast between all those calculated curves of different compositions. p(T) curves were 
obtained by fitting experimental data with Eq. (32). Table 1 lists all the fitting parameters namely, T-independcnt 
scattering rate constant (A2), charge carriers concentrations (y/np), charge gap parameter (Ei) and dp(T)/dT in 
detail above T c and T*. ^Jnp is estimated by first computing Ej from Eq. (32) and then substitute it into Eq. (28). 
Furthermore, two values of ^Jnp are computed at 300 K, in which one satisfies m* /m = 1 while the other takes m* /m 
= 50. m is the rest mass of an electron and the fitted values of A 2 are actually equal to Aziih 2 / e 2 k,B- Sulkowski et 
al. [31], have computed charge carriers concentrations from Hall effect measurements for TmE^CuaOs.gg that is in 
the order of 10 22 cm~ 3 . This value is comparable with the estimated value from Eq. (28) with m* /m — 50, which is 
4.7 x 10 18 m~ 2 = 4.7 x 10 22 cm~ 3 . Parallel to this, Eqs. (32) and (28) also enable one to approximate the effective 
mass of charge carriers. As for single crystals, one may still employ Eq. (32) via the assumption [16] 



P P oiy(T) = Vpc(T)p ab (T). (43) 

Eq. (43) is assumed to be valid to convert both p c (T) and p a b(T) data to p po i y (T), which will enable one to 
analyze the effect of A2 parameter between single crystals and polycrystals. One should not substitute Eqs. (33) 
and (34) into Eq. (43) because Eq. (34) was derived based on a special case as mentioned previously. Subsequently, 
the changes on scattering rate (A2) will be captured. Recently, Lanzara et al. [32] have shown quite convincingly via 
ARPES measurements that the existence of electron-phonon coupling associated with movements of oxygen atoms in 
Bi 2 Sr 2 CaCu20 8 , Bi 2 Sr 2 Cu20 6 and La 2 _ a; Sr :l; Cu04 should not be neglected entirely. This observation could be due 
to polarons that is well represented by Ej in iFDS. The important difference between polarons and free electron- 
phonon scattering is that the latter has a very strong T-depcndence while the former increases the effective mass of 
the charge carriers to some extent. This could be the sole reason for the missing electron-phonon coupling effect on 
p(T) measurements in high-T c superconducting cuprates thus far. All Ej values were obtained from Ref. [33] and 
the predictions stated above are only valid for reasonably pure materials without any significant impurity phases and 
grain boundary effects. 

V. CONCLUSIONS 

In summary, iFDS is shown to be an appropriate theory in the normal state of high-T c superconductors. The varia- 
tions on polycrystalline resistivity with doping were very well explained with one exception, in which the importance 
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FIG. 1: Computed p(T, Ei) curves for Hgo.7Cro.3Sr2Cu04+5 and Gdo.95Pro.osBa2Cu307-i from Refs. [29] and [26] respectively. 
The inset is for Sro.gLao.iCuO^ and TmBa2Cu306.99 from Refs. [30] and [31] respectively. The curves were determined in 
accordance with Eq. (32). The symbols were used solely to distinguish curves of different samples and it does not represent 
experimental data points. Apparently, T cro3SO „ e rS for all these samples (if any) are < T c s, as such the metal to semiconductor 
transition is not observed, which has been anticipated because Ei < ksTc. 



TABLE I: Computed values of the temperature independent scattering rate constant (A2), charge gap parameter (Ei), resistivity 
slope (dp(T)/dT) and the concentration of charge carriers (y / np) at 300 K with different effective masses (m*/m = 1 and 
m*/m = 50) in the normal state of high-T c superconducting polycrystals. Note that the magnitudes of Ei are in the range of 
2-9 meV and the charge carriers' concentration is determined to be in the order of 10 18 m~ 2 or 10 22 cm~ 3 if m*/m = 50 is 
assumed. dp(T)/dT is calculated from the experimental p po i v {T) linear plots. 

of characteristics temperature was ignored throughout this work partly due to vague understanding on temperature 
dependence of c-axis and a6-plane scattering rates. It is interesting to note that iFDS is able to predict the resistivity's 
magnitude with doping via ionization energy as well as enable one to estimate the valence state of substituting ions 
accurately. Apart from that, the fitting parameters give a reasonable prediction of the charge carriers' effective mass in 
the vicinity of 10 -29 kg or 50 times heavier than electron's rest mass. The concentration of this heavy charge carriers 
is around 10 22 cm~ 3 whereas light charge carriers' concentration is approximately 10 20 cm -3 when the effective mass 
is equals to rest mass. The heavier effective mass (1 < m* /m < 50) somewhat imply the existence of polarons in 
the normal state of high-T c superconductors that could mask electron-phonon coupling effect on resistivity versus 
temperature measurements. 
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Samples, [Ref.] 


A 2 


E l 
[K] 


Ei 
[meV] 


T^pxlO 16 

[m -2 ] 300 K 
m*/m - 1 


[m 2 ] 300 K 
m*/m - 50 


dp(7) / dJ 
[liQ-cmK 1 ] 


TmBa 2 Cu 3 6 .99 [31] 


1.6455 xlO" 3 


43.1 


3.71 


9.354 


4.677 


1.60 


Sr .9Lao.iCu0 2 [30] 


9.6000 x 10 4 


47.8 


4.20 


9.208 


4.604 




Hg .7Cro.3Sr 2 Cu0 4+ 5 [29] 


1.0250 x Iff 2 


30.0 


2.59 


9.771 


4.886 


1.00 


Gdo.95Pro.o5Ba 2 Cu30 7 _s [26] 


3.8200 x 10" 3 


94.4 


8.13 


7.884 


3.942 


3.10 


YBa 2 Cu 3 7 Sample B, [16] 


1.1458x10* 


89.2 


7.69 


8.022 


4.011 




YBa 2 Cu 3 7 Sample C, [16] 


1.1268 x Iff 8 


94.0 


8.11 


7.895 


3.948 




GaSr 2 (Yo. 6 Ca .4)Cu 2 7 , [16] 


1.3002 xlO" 5 


69.3 


5.98 


8.572 


4.286 


~0 


YSrBaCusO? Sample [O], [16] 


1.1712xl0 6 


53.5 


4.61 


9.036 


4.518 


1.10 


Y .8Pro. 2 SrBaCu 3 7 Sample [0], [16] 


1.8589 xlO" 6 


76.6 


6.61 


8.366 


4.183 


1.60 


YBa 2 Cu 3 O x , [16] 


5.5747 x 10" 6 


95.6 


8.25 


7.853 


3.926 


4.25 


YBa 2 Cu 3 7 . x Sample A, [16] 


6.3767 x 10" 6 


92.0 


7.94 


7.948 


3.974 


5.00 


Y(Ba .85Pbo.i5)2Cu 3 7 _ x Sample B, [16] 


5.4926 x 10" 6 


79.8 


6.88 


8.278 


4.139 


4.90 


Tl(Sr 1 . 6 Sm .4)CaCu 2 O 7 . x , [16] 


2.8695 x 10" 5 


59.5 


5.13 


8.857 


4.430 


2.60 


Tl(Sr 1 . 6 Dy .4)CaCu 2 O 7 . x , [16] 


2.2192 x 10" 5 


66.2 


5.71 


8.661 


4.330 


2.00 



